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Abstract 

One of the difficulties encountered when studying physical theories 
in discrete space-time is that of describing the underlying continuous 
symmetries (like Lorentz, or Galilei invariance). One of the ways of 
addressing this difficulty is to consider point transformations acting 
simultaneously on difference equations and lattices. In a previous 
article we have classified ordinary difference schemes invariant under 
Lie groups of point transformations. The present article is devoted to 
an invariant Lagrangian formalism for scalar single- variable difference 
schemes. The formalism is used to obtain first integrals and explicit 
exact solutions of the schemes. Equations invariant under two- and 
three- dimensional groups of Lagrangian symmetries are considered. 



1 Introduction 



A recent article was devoted to a symmetry classification of second order ordi- 
nary difference equations [1]. This was modeled on a paper by Sophus Lie, in 
which he provided a symmetry classification of second order differential equations 
(ODEs) [2]. As a matter of fact, the classification of difference schemes goes over 
into Lie's classification of ODEs in the continuous limit P]. 

S. Lie showed that a second order ODE can be invariant under a group Gr 
of dimension N = 0,1, 2, 3, or 8. For N > 2 the equation can be integrated 
in quadratures. This can be done by transforming the equation to one of the 
"canonical" forms, integrated by Lie himself [2]. Virtually all standard methods 
of integrating second order ODEs analytically can be interpreted in this manner 
(though this is not mentioned in most elementary textbooks). 

The situation with difference equations is much less developed. This is not 
surprising, since applications of Lie group theory to difference equations are much 
more recent [3,. ..,27]. Several different approaches are being pursued. One possi- 
bility is to consider the difference equations on a fixed lattice [3,. ..,13] and consider 
only transformations that do not act on the lattice. In order to obtain physically 
interesting symmetries in this approach, it is necessary to go beyond point sym- 
metries and to let the transformations act on more than one point of the lattice. 
Lie algebra contractions occur in the continuous limit and some "generalized" 
symmetries may "contract" to point ones |lUj . 

The second possibility is to consider group transformations acting both on 
the difference equations and on the lattice [1,17,. ..,27]. Technically, for systems 
involving one dependent and one independent variable, this is achieved by con- 
sidering a difference scheme, consisting of two equations, one representing the 
actual difference equation, the other the lattice. 

This is the approach that we will follow in the present article. More specifi- 
cally, we will consider the same three-point scheme as in our previous article PP. 
The continuous limit of the scheme, if it exists, will be a second order ODE. 

Thus, we consider two variables, x and y, with x the independent one and y 
dependent. The variable x runs through an infinite set of values {x = Xk, k E 1^} 
that are not necessarily equally spaced and are not prescribed a priori. Instead, 
we give two relations between any three neighboring points 



and also specify some initial conditions like Xq, Xi, yo = yi^Xo), yi = y{xi). In the 
continuous limit eq. ()1.1|) goes into an ODE, ()1.2j) into an identity (like = 0), 
if the continuous limit exists. 
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The group transformations considered in this approach are of the same type 
as for ODEs. They are generated by a Lie algebra of vector fields of the form 



X = ax,y)^ + Vix,y)^. (1.3) 

The corresponding transformations are purely point ones, since the coefficients ^ 
and 7] depend only on {x,y), not on the shifted points or 

In we showed how Lie group theory can be used to classify such pairs 
of equations as and ()1.2j) . Possible dimensions of the symmetry group G 
of eq. jHH), are N = 0, 1, 2, 3, 4, 5, and 6. The highest dimension, = 6, 
occurs only for difference schemes equivalent to 

— ^ \^ =0, x+ — 2x + X- = 0. 

(x+ — x)^ 

The purpose of this article is to provide a Lagrange formalism and difference 
analog of Noether's theorem for second order difference schemes of the form p.ll) 
and p.2|l . admitting Lie point symmetry groups. The Lagrangians will be used 
to obtain first integrals and exact analytic solutions of the difference schemes. 



2 General theory 

2.1 Definitions and notations 

We study the difference system (jl.lj) and ()1.2|1 . In general, we assume that these 
equations can be solved to express x+ and y+ explicitly in terms of (x, y, X-,y-) 
and also vice verse, i.e. in terms of the other quantities. We also make 

use of the following quantities 

, _ , _ _ y+-y y-y- 

~ x+ ~ X, — X — X-, yx — — , y: 



2 ^^'^^ 

yxx ~i ; ; ijjx 
+ h_ 

i.e. the up and down spacings in x, the right and left discrete first derivatives 
and the discrete second derivative, respectively. It is also convenient to use the 
following total shift and discrete differentiation operators: 

S - 1 

5/(x) = /(x±), D = ^!^. 
±h ±h ±h± 

Continuous first and second derivatives are denoted y' and y", respectively. 
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When acting on differential equations, the vector fields ()1.3p must be pro- 
longed to act on derivatives. For difference schemes, the prolongation of a vector 
field acts on variables at other points of the lattice. It is obtained by shifting the 
coefficients to the corresponding points. For three point schemes we have 

d d d d 

prX = X + ^{x_,y_)- — + ^{x+,y+)- — +r]{x-,y-)^— + r]{x+,y- ^ 



ox^ ox^ oy- oyj^ 

(2.2) 



2.2 Lagrangian formulation for second order ODEs 

It has been known since E. Noether's fundamental work that conservation laws 
for differential equations are connected with their symmetry properties [28,..., 31]. 
For convenience we present here some well-known results adapted to the case of 
second order ODEs. 

Let us consider the functional 

Uy) = J^Lix,y,y')dx, IcR\ (2.3) 

where L{x, y, y') is called a first order Lagrangian. The functional ()2.3|) achieves 
its extremal values when y{x) satisfies the Euler-Lagrange equation 

6L _ dL _ ^ rdL\ _ ^ d ,d „d 

6y dy \dy' ) ' dx ^ dy ^ dy' 

where D is the total derivative operator. The equation ()2.4|1 is an ODE that can 
be rewritten as 

y" = f\x,y,y'). (2.5) 

Let us consider a Lie point transformation G generated by the vector field 
()1.3|1 . The group G is a "variational symmetry" of the functional L if and only 
if the Lagrangian satisfies 

prX(L) + LD{C) = 0, (2.6) 

when prX is the first prolongation of the vector field X for y'. We will actually 
need a weaker invariance condition than given by eq. ()2.(j|l . The vector field X 
is an "infinitesimal divergence symmetry" of the functional C{y) if there exists a 
function V{x, y) such that [28] 

wX{L) + LD{i) = D{y), V = V{x,y). (2.7) 



The two important statements for us are: 

1. If X is an infinitesimal divergence symmetry of the functional £, it gen- 
erates a symmetry group of the corresponding Euler-Lagrange equation. The 
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symmetry group of eq. ()2.4|) can of course be larger than the one generated by 
symmetries of the Lagrangian. 

2. Noether's theorem [28,.. .,31] can be based on the following Noether-type 
identity [31], which holds for any vector field and any function L: 

prX(L) + LD{i) = (r^ - iy')— + D{iL + (r? - iv')^)- (2-8) 

It follows that if X is a divergence symmetry of L, i.e. ()2.6p or ()2.7|) is satisfied, 
then there exists a first integral 

BT 

^L + {r]- iy')— -V = K = const (2.9) 
dy' 

of the corresponding Euler-Lagrange equation. 

The above considerations tell us how to obtain invariant ODEs and conser- 
vation laws from divergence invariant Lagrangians. They do not tell us how to 
obtain invariant Lagrangians for invariant equations. This amounts to "varia- 
tional integration" , as opposed to variational differentiation. 

A procedure that we shall use below to find invariant Lagrangians for differential 
equations can be summed up as follows. 

Start from a given ODE y" = f{x,y,y') and its symmetry algebra with basis 

d d 
Xa = ^a{x,y)— + T]a{x,y)—, a = l,...,k. 
ox ay 

Find the invariants of Xa in the space {x,y,y',A}, where A is the Lagrangian. 
The appropriate prolongation in this case is 

prX = ^^ + V^y+e^- (DOA-^, C = Di,) - y'DiO (2-10) 

and we require that L{x,y,y') should satisfy 

prX(A-L)|A=L = 0. (2.11) 
Each basis element X^ provides us with an equation of the form 

Solve the partial differential equations ()2.12|) . This will give us the general 
form of an invariant Lagrangian. It may involve arbitrary functions of the invari- 
ants of X. 

Request that the Euler-Lagrange equation (j2.4|) should coincide with the equa- 
tion we started from. This will further restrict the invariant Lagrangian and 
determine whether one exists. 
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If this procedure does not yield a suitable Lagrangian, then step 1 can be 
weakened. We can request that the Lagrangian be invariant under some subgroup 
of the symmetry group of the given ODE, rather then the entire group. We then 
go through step 2, then verify whether the obtained Lagrangian is divergence 
invariant under the entire group, or at least a larger subgroup. In any case, each 
divergence symmetry of the Lagrangian will provide a first integral of the ODE. 

For ODEs the Lagrangian formalism is not the only integration method. The 
existence of one-parameter symmetry group provides a reduction to a first-order 
ODE directly. The existence of a two-parameter symmetry group makes it possi- 
ble to integrate in quadratures. An invariant Lagrangian provides an alternative. 
Indeed, assume that we know two first integrals 

h{x,y,y') = A, Mx,y,y')=B, (2.13) 

then we eliminate y' from these two equations and obtain the general solution 

y = F{x,A,B), (2.14) 

of the corresponding ODE ()2.5j) by purely algebraic manipulations. It is this 
method of invariant Lagrangians that generalizes to difference equations and is 
particularly useful when direct methods fail. 

2.3 Lagrangian formalism for second order difference equa- 
tions 

The variational formulation of discrete equations and a discrete analog of Noether's 
theorem are much more recent [19,25,26,27]. Here we briefly overview the results 
that we shall need below. 

Let us consider a finite difference functional 

hh = Y^ C{x, x+, y, y+)h+, (2.15) 

defined on some one-dimensional lattice f2 with step /i+ that generally can depend 
on the solution 

h+ = ^{x,y,x+,y+). (2.16) 

The Lagrangian (|2.15|) must be considered together with a lattice (2.16). On 
different lattices it can have different continuous limits and in this limit the lattice 
equation itself vanishes (turns into an identity like = 0) 

In the continuous case, a Lagrangian L provides an equation (the Euler- 
Lagrange equation) that inherits all the symmetries of L. In the discrete case we 
wish the Lagrangian (j2.15|) to provide two equations: the entire difference sys- 
tem (1.1), (1.2). Moreover, the three-point difference system should inherit the 
symmetries of two-point Lagrangian. 
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Let us again consider a Lie group of point transformations, generated by a 
Lie algebra of vector fields of the form p.3|) . The infinitesimal invariance 
condition of the functional ()2.15jl on the lattice ()2.16j) is given by two equations 
[19,25,27]: 

^ (2.17) 



-h 



where 



e+=e(a:+,y+), v^ = v{x+,y+)- (2.18) 

Let us consider a variation of the difference functional ()2.15j) along some curve 
y = (pix) at some point {x,y). The variation will effect only two terms in the 
sum fjTlb : 

L,, = ... + C{x, x_,y, y-)h_ + C{x, x+, y, y+)h+ + (2.19) 
so we get the following expression for the variation of the difference functional 

5L = + ^-fsy, (2.20) 
dx by 

where 5y = 4>'6x and 

5L BC , dC- ^_ ^ 5C , dC , dC' 

— = h+— + h^— +C - C, — = h+— + h^— , 
ox ox ox oy Oy Oy 

where C~ = S (C). 

Thus, for an arbitrary curve the stationary value of difference functional is 
given by any solution of the two equations, called quasiextremal equations 

ox oy 

Both of them tend to the differential Euler-Lagrange equation in the continuous 
limit. Together they represent the entire difference scheme and could be called 
"the discrete Euler-Lagrange system". The difference between these two equa- 
tions, or some other function of them that vanishes in the continuous limit will 
represent the lattice. 

Now let us consider a vector field ()1.3|) with given coefficients ^{x,y) and 
ri{x,y). Variations along the integral curves of this vector field are given by 
6x = ^da and 6y = rjda, where da is a variation of a group parameter. A 
stationary value of the difference functional ()2.15|1 along the flow generated by 
this vector field is given by the equation: 
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which depends exphcitly on the coefficients of the generator . 

If we have a Lie algebra of vector fields of dimension 2 or more, then a station- 
ary value of the difference functional ()2.15|1 along the entire flow will be achieved 
on the intersection of the solutions of all equations of the type ()2.22|1 . i.e. on the 
quasiextremals (j22H)- 

On the other hand, equations ()2.21|) can be interpreted as a three-point dif- 
ference scheme of the form fll.l|) . ()1.2|) . For instance, given two points {x,y) and 
(a;_,?/_), we can calculate In the continuous limit both of these equa- 

tions will provide the same second-order differential equation. Thus, one of the 
quasiextremal equations can be identified with eq. (jl.lj) and the difference be- 
tween the two of them with the lattice equation ()1.2|) . 

It has been shown elsewhere [19,25,27], that if the functional ()2.15|) is invariant 
under some group G, then the quasiextremal equations 1)2.211) are also invariant 
with respect to G. As in the continuous case, the quasiextremal equations can 
be invariant with respect to a larger group than the corresponding Lagrangian. 

A useful operator identity, valid for any Lagrangian C{x,x^,y,y^) and any 
vector field X is ([19,25]): 

/ dC- , dC- 
+D h^T]— +h.^— 
+h \ oy ox 

From eq. (|2.23|) we obtain the following discrete analog of Noether's theorem. 

Theorem 2.1 Let the Lagrangian density C be divergence invariant under a Lie 
group G of local point transformations generated by vector fields X of the form 
()1.3|) . i.e. let us have 

prXiC) + CDiO = DiV) (2.24) 

+h +h 

for some function V{x, y). Then each element X of the Lie algebra corresponding 
to G provides us with a ffist integral of the quasiextremal equations ()2.2ip . namely 

K = h^r]— +h^^— +^C--V. (2.25) 
oy ox 

Proof. [19,25] On solutions of the quasiextremal equations ()2.2ip eq. ()2.23p re- 
duces to 

i,(Mf^Mf^?£-)^5W (2.26) 

(we have used eq. (|2.23|) ). The result (|2.25p follows immediately. □ 
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The fundamental equation ()2.25p is the discrete analog of eq. ()2.9p for ODEs. 

Let us compare the situation for second order ODEs and for three-point differ- 
ence schemes. For a second order ODE a Lagrangian that is divergence invariant 
under a two-dimensional symmetry group provides two integrals of motion. From 
them we can eliminate the remaining first derivative and obtain the general so- 
lution, depending on two arbitrary constants (the two first integrals). Moreover, 
we do not really need a Lagrangian. Once we have a two dimensional symmetry 
group of the ODE, we can integrate in quadratures. 

For three-point difference schemes we have two equations to solve, namely the 
system (1.1), (1.2). Equivalently, we have a set of points (xn,?/„), labeled by an 
integer n. Any 3 neighboring points are related by two equations that we can 
write e.g. as 

yn+l = Fi{Xn,yn,Xn-l,yn-l), Xn+1 = ^l{Xn, Vn, Xn-1, Vn-l) ■ {2.27) 

Alternatively, the system could be solved for i/n-i- We mention that we use 
notations like x„_i = X-, Xn = x, Xn+i = x+, yn-i = y-, i/n = y, yn+i = y+ 
interchangeably. 

Given some starting values {xo,yo,X-i,y-i), we can solve ()2.27|1 for (x„,y„) 
with n > 1, and n < —2. The solution will depend on four constants Ki,i = 
1, ...,4, and can be written as 

yn = yn{xn, i^i, K2, K3, K^), (2.28) 

Xn = Xn{Ki,K2,K3,K^). (2.29) 

The two quasiextremal equations ()2.21|) correspond to the system fj2.27p . 

A one-parameter symmetry group of the Lagrangian C will provide us with a 
first integral (|2.25|) . i.e. an equation of the form 

f{Xn, Vn, Xn+1, l/„+l) = Ki. (2.30) 

compatible with the system ()2.21|) . We can solve 1)2.301) for e.g. yn+i, substitute 
into ()2.27|) and thus simplify this system. 

A two-dimensional symmetry group will provide two first integrals of the 
form ()2.25j) . We can solve for Xn+i and yn+i- Then system (I2.27|l is reduced to a 
two-point difference scheme. Quite often it is possible to solve it by integration 
methods that allow one to integrate a two-point difference scheme explicitly. 

A three-dimensional symmetry group provides three first integrals of the type 
fl2.25|) . From them we can express Xn-i, Vn-i and yn in terms of x„. This provides 
us with the solution ()2.28|) and a two-point difference equation relating Xn+i and 
Xn- If this equation can be solved, we have a complete solution of the problem. 
Finally, if we have four first integrals, then we get the general solution of the 
system by purely algebraic manipulations. 
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An alternative method can be proposed when the Lagrangian is invariant with 
respect to a two-dimensional Lie group. The discrete Lagrangian corresponding 
to a given continuous one is not unique and it is possible to introduce a family 
of Lagrangians: 

Ci = Ci{x,x+,y,y+,ai,Pi), i = 1,2,3,... (2.31) 
depending on parameters ai. Pi, all satisfying 

lim(^x+,y+)-*{x,y)^i{x, x+, y, y+, ai, pi) = C{x, y, y') 

for the same continuous Lagrangian C{x,y,y'). 

Let us take three different Lagrangians in the family ()2.3H1 . corresponding to 
constants ai, Pi,a2, P2 and a^, P^. Each of them will lead to two first integrals 
and two quasiextremals. In examples considered below we will show that it is 
possible to fine-tune the constants ai,Pi in such a manner as to get a system of 
two invariant equations of the form ()l.H) . ()1.2j) and three first integrals, yielding 
a set of solutions to the two quasiextremal equations. It is these two equations 
that will constitute the invariant difference system. 

In Section 2.2 we described a procedure for obtaining invariant Lagrangians 
for given second-order differential equation. For difference equations our starting 
point will be a discretization of the continuous Lagrangian. This is obviously not 
unique and we shall make use of the inherent arbitrariness. Once an invariant 
difference Lagrangian with a correct continuous limit is chosen we construct the 
invariant difference scheme in the manner described above. 

In our previous article [1] we gave a classification of difference schemes and 
used all realizations of Lie algebras that provide such schemes. Any algebra 
containing a two-dimensional subalgebra realized by linearly connected vector 
fields such as 

f d_ d_\ fd_ d_\ 
\dx^ ^ (9x / ' \ (9x ' (9x / ' 

leads to a linear differential equation and its discretization. 

Below we shall consider only genuinely nonlinear difference schemes presented 
in Ref P that have nonlinear differential equations as their limit. 

3 Equations corresponding to Lagrangians invariant un- 
der one and two-dimensional groups 

3.1 One-dimensional symmetry group 

We start with the simplest case of a symmetry group, namely a one-dimensional 
group. Its Lie algebra is generated by one vector field of the form ()1.3|1 . By an 
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appropriate change of variables we take this vector field into its rectified form. 
Thus we have 

D,, : X. = |. (3.1) 

The most general second order ODE invariant under Xi is 

y" = F{x,y'), (3.2) 

where F is an arbitrary given function. 

Eq. ()3.2|) is actually already reduced to a first order equation for u = y' . If 
Xi is a variational symmetry of eq. (j3.2j) and we know the Lagrangian L that it 
comes from, we can do better. An invariant Lagrangian density will by necessity 
have the form L = L{x,y') (see eq. ()2.6|) ). The Euler-Lagrange equation ()2.4p 
reduces to 

Substituting for y" from eq. ()3.2j) . we obtain a linear partial differential equa- 
tion for L{x,y'). This of course has an infinity of solutions. Let us assume that 
we know a solution L{x,y) explicitly. Eq. ()2.9|) . i.e. Noether's theorem, provides 
us with a first integral 

8T 

Q^,M)=K. (3.4) 

We can solve (in principle) eq. ()3.4j) for y' as a function of x (and K). The general 
solution is then obtained by a quadrature: 

y' = ^{x,K), y(x) = yo+ <P{x,K)dt. (3.5) 



In the discrete case the situation is similar. Let us assume that we know 
a Lagrangian C{x,x+,y,y+), invariant under the group of transformations of y, 
generated by Xi of eq. ()3.1|1 . It will have the form 

C = C{x, x+, y^), y^ = — — -. (3.6) 

x+ — X 

The corresponding quasiextremal equations, to be identified with the system (jl.lj) 
and p.2j) . are 

5C dC , , dC , \ ^ /o ^,\ 

^ =-^(^' y^) + — (^-' y^) = 0; (3-7) 

dy oy^ yx 

6C _ dC dC \ ri \ 

'T~ — '^+T^(3;, a;_|_, yx) + yx-^ (2;, yx) — L-yx, a;_|_, yx) 
ox ox oyx 

(3.8) 

+h_ — [x_,x,yx) - yx-g;i^{x-,x,yx) + £(x_, x, y^.) = 0. 
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The first integral ()2.25|) can be read off form eq. ()3.7|) and is 

dC 

— {x,x+,y^) = K. (3.9) 
dyx 

We can solve eq. ()3.9|) for and by down shifting obtain y^: 

ycc = 4>{x,x+,K), ys: = 4>{x-,x,K). (3.10) 

Substituting into the quasiextremal equation ()3.8|) . we obtain a relation between 
x^, X- and x, i.e. a single three-point relation for the variable x. For y we then 
obtain a two point equation 

y+-y = {x+-x)(f){x,x+,K). (3.11) 

Equation (j3.1ip is really a discrete quadrature: a first order inhomogeneous linear 
equation for y. 

Example 3.1 Consider the Lagrangian 

£ = <x^+,exp(y,.). (3.12) 
The quasiextremal equations are 

exp(?/^.) + exp(?/5.) = 0; (3.13) 



h+ax'^ ^x^+i exp(?/^) + y^x^x''^^^ exp{y^) - xlx\^^ exp(?/^) 

exp(?/s) - i/s<_ia;^ exp(ys) + x^.^x^ exp(ys) = 0. 

The first integral is 



(3.14) 



x>^+iexp(y,) = K (3.15) 
From ()3.15|) we have 

y, = HKx-^x-l), y, = \n{Kx-l,x-'). (3.16) 

Eq. ()3.13|) is satisfied identically. Eq. ()3.14|) reduces to a three point equation for 
x: 

ax„+i + (6-a)x„-te„_i+x„(-61n(x„+i) + (6-a) ln(a;„)+aln(x„_i)) = 0. (3.17) 
This lattice equation can be reduced to a two-point equation for a new variable 

b 

aXn + {b-a) - = aln(A„_i) + 61n(A„). (3.18) 

12 



In particular, one can choose the solution = A„+i = A, where A satisfies the 
equation 

h 

aX + {b-a) - - = (a + b) ln(A). 
A 

It provides us with the lattice x„ = xqA". 

Substituting the lattice into ()3.15p . we obtain a two point equation for y: 

Vn+i -yn = (xn+i - x^) ln{K x""" x"^^) . (3.19) 

The fact that we could solve eq. ()3.17p explicitly is specific for the considered 
example. The fact that we obtain a three-point equation involving only the 
independent variables is true in general. 



3.2 Two-dimensional symmetry groups 

J^2,i The Abelian Lie algebra with non-connected basis elements 

corresponds to the invariant ODE 

y" = Fiy'), (3.21) 

where F is an arbitrary function. 

The equation can be obtained from the Lagrangian 

L = y + G{y'), Hv') = (3-22) 

The Lagrangian admits symmetries Xi and X2: 

piXiL + LD{^i) = 0; 

prX2L + LD(e2) = l = D{x). 

With the help of Noether's theorem we obtain the following first integrals: 

Ji = y + G{y') - y'G'iy'), = G'{y') - x. (3.23) 

As we mentioned in the Section 2.2, it is sufficient to have two first integrals to 
write out the general solution of a second order ODE without quadratures. More 
explicitly, we can solve the second equation ()3.23|) for y' in terms of x and obtain 

y' = H{J2 + x), H{J2 + x) = [G']-\J2 + x). (3.24) 
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Substituting into the first equation, we obtain 

y{x) = Ji- G[H{J2 + x)] + (J2 + x)H{J2 + x). 



(3.25) 



Now we are in a position to show how one can find a variational discrete model 
and its conservation laws by means of Lagrange-type technique. Let us choose a 
difference Lagrangian in the form 

C = y^ + G{y^), (3.26) 

then 

pvXiC + CD{ii) = 0; 

+h 

PVX2C + CDi^2) = 1 = Dix). 

+h +h 

The variations of C yield the following quasiextremal equations: 

^: G'{y.)-G'{y,) = ^^±±^; (3.27) 

Due to the invariance of the Lagrangian with respect to the operators Xi and 
X2, the difference analog of Noether's theorem yields two first integrals 

h = y + G{y^) - y.,G\y,) + ^^^i/., (3.29) 

h = G\y.^) - (3.30) 
As in the case of the algebra Di^i we can solve for y^ to obtain 

= $i(J2,x + x+). (3.31) 
Substituting into the equation for Ji we obtain 

7/ = <l>2(/l,/2,X,X+). (3.32) 

Calculating y^ from eq. ()3.32p and setting it equal to ()3.3H) . we obtain a three 
point recursion relation for x. Solving it (if we can), we turn eq. ()3.32j) into an 
explicit general solution of the difference scheme ()3.27p . ()3.28p . 

Example 3.2 Let us consider the case 

^ = ^"^2^ + exp(?/^). (3.33) 
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The two first integrals ()3.29|) . ()3.30|) in this case are the following: 

Ii = y + exp(?/^.) - exp(?/^) + ^"^ y^^ (3.34) 



h = exp(?/^) - 



2 

Equation ()3.31|) and ()3.32|) reduce to 



y. = ln[h + ^^^^^], (3.35) 



y = h-h- ^^^^±1±^ + ih + x„) In i^h + ^!^±1±^ j , (3.36) 
The recursion relation for x is 

^"'^^^"^ + {h + X„)[ln(2/2 + Xn+l + Xn) - ln(2/2 + Xn+ Xn-l)]- (3.37) 

The last equation is difficult to solve. We have however reduced a system of two 
three-point equations to a single three-point one. We shall return to this case in 
Section 5 using an alternative method. 

D2,2 The non-Abelian Lie algebra with non-connected elements 

d d d 

Xi = — , X, = x— + y— (3.38) 
ay ox ay 

yields the invariant ODE 

y"=-F{y'). (3.39) 
We define a function G{y') by the equation 

F(y-) - (3.40) 

Then the ODE ()3.39|1 is the Euler-Lagrangian equation for the Lagrangian 

L = -G{y'), (3.41) 

X 

which admits Xi and X2 as variational symmetries: 

piXiL + LDi^i) =0; 

PtX2L + LD{^2) =0. 
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Noether's theorem provides us with two first integrals: 



Ji = -G'(y'), J2 = G{y') + - y') G\y'). 
X \x / 

Let us take the difference Lagrangian 

X ~r X-i- 

which satisfies 

prXiC + CDiCi) =0; 

+h 

prX2£ + £D(6) =0. 
Then the variations of C yield the following quasiextremal equations: 

oy X + x+ X- + X 

2h 2 2 

^ ^ .^ Giy.) + T- — ^G'iy.)yx - ^— — ,Giy.) (3.42) 

^£ + [x + x+j [x + x+j 

' 2/i_ . 2 . 2 



- . —G'{y,)y, + —My,) = 0. 



(x_ + xy (x_ + x) (x_ + x) 

Since the Lagrangian is invariant with respect to the operators Xi and X2, 
we find the first integrals 

2G'(yJ , Axx. ^, , 2G'(y^ , , , , , 

X ~r IX ~r x^ i x ~r 

for the solutions of ()3.42|) . 

As in the case of the algebra D2,i, we can solve for y^, using the integral Ii. 
We obtain: 

y^ = ^i{Ii,x + x+). (3.44) 
The second integral allows us to express y as a function of x and x+ 

^ = x<|.i + ^ - G($0. (3.45) 
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4 Equations corresponding to Lagrangians invariant un- 
der three-dimensional Lie groups 

Among the "prototype equations" of our previous article [1], many have three- 
dimensional symmetry groups. In this section we shall consider two of these cases. 
Both of them come from Lagrangians that also have three-dimensional symmetry 
groups, i.e. all these symmetries are Lagrangian ones. 

D3 1 Let us first consider a family of solvable Lie algebras depending on one 
constant k: 

Xi = — , ^2 = — , X3 = x— + %— , A;^0,-,l,2. (4.1) 
ox ay ox ay 2 

The invariant equation has the form 

y" = y'^^. (4.2) 

This equation can be obtained by the usual variational procedure from the 
Lagrangian 

L=^-^{y')—i+y, (4.3) 
which admits operators Xi and X2 for any parameter k and X3 for = — 1: 

prXiL + LD{^i) = 0; 

PtX2L + LD{^2) = 1 = D{x)- 

pTX,L + LD{^,) = {k + l)L. 

It is possible to show that there is no Lagrangian function L{x,y,y') which 
gives eq. ()4.2|) with A; 7^ —1 as its Euler's equation and is divergence invariant for 
all three symmetries ()4.1|) . 

For arbitrary k there are two first integrals 

Ji = +y = A', J2 = {k- l)(y')^ -x = B\ 

Eliminating y' we find the general solution: 

In the case /c = — 1 we have the further first integral corresponding to the 
symmetry X3: 

h = ^{y-xy')+xy = C''. 
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It is functionally dependent on Ji and J2 since a second order ODE can possess 
only two functionally independent first integrals. Let us mention that the first 
integral J3 is basic: 

Jl=Xi(J3), J2 = -X2(J3), 

since 

[Xi,X3]=Xi, [X2,Xs] = kX2. 

In this case we have the following relation: 

4-JiJ2-J3 = 0. (4.5) 

Thus, the integral J3 is not independent and is of no use in the present context. 

Now let us turn to the discrete case and consider k = — 1 only. Other values 
of k will be considered in Section 5, using a different approach. Let us choose the 
Lagrangian to be 

C = -A^,+ y^ (4.6) 

as a discrete Lagrangian, which is invariant for Xi and X3 and divergence invari- 
ant for X2: 

piXiC + CDi^i) = 0; 

+h 

+ CD{i2) = l = D{x)- (4 7) 

+h +h ^ > 

+h 

From the Lagrangian we obtain the quasiextremal equations: 
5C 4/1 1 



1: 



(4.^ 



w ^ y + y+.y- + y n 



This system of equations is invariant with respect to all three operators ()4.H1 . 
The application of the difference analog of the Noether theorem gives us three 
first integrals: 

J o ^ , y + y+ A T 2 X + X+ 

h = -2^ + — - — = A, l2 = — — — = B, 

2{x+y-y+x) x+y + y+x ^ ' > 

-'3 - 1 'F= ^ n ~ 

In contrast to the continuous case the three difference first integrals Ii , I2 and 
/3 are functionally independent and instead of eq. ()4.5j) we have the following 
relation: 

4-^/2-/3 = ^/^^?/.= (^f^. (4.10) 
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This coincides with eq. ()4.5|) in the continuous hmit e 0. We see that the ex- 
pression h'^y^ is also a first integral of ()4.8|) . This allows to introduce a convenient 
lattice, namely: 

1 1 4e^ 

-^h'iyx = -^h\y^ = (^^^2y ' ^ " < £ < 1. (4.11) 

Substituting from eq. 1)4.111) into I2, we obtain a two-term recursion relation 
for X, namely 

Xn+i-{l + e)xn-eB = 0, (4.12) 

or 

-{1 + e)xn+i + Xn - eB = 0, (4.13) 

depending on the sign choice for ^/y^. These equations can be solved and we 
obtain a lattice satisfying 

x„ = (xo + 5)(l+£)"-5, xo>-B (4.14) 

for the first equation and a lattice satisfying 

Xn = {xo + B){l + e)-" -B, xo<-B (4.15) 

for the second equation. Using the expressions for Ji, we get the general solution 
for y (the same for both lattices ()4.14j) and ()4.15j) ) as 

This agrees with the continuous case up to order e. 

We have used the three integrals Ji, I2 and I3 to obtain the general solution of 
the difference scheme ()4.8|) . Indeed, the solution ()4.14|) . ()4.16|) for x„, yn depends 
on 4 constants {A, B,xo,e), as it should. 

The difference scheme is not compatible with a regular lattice, but requires 
an exponential one, as in eq. ()4.14j) . The only non-algebraic step in the inte- 
gration was the solution of eq. ()4.12j) : a linear two point equation with constant 
coefficients. 

D3,2 The group given by the operators 

Xi = — , X2 = 2x— + y—, Xs = x'—+xy— (4.17) 
ox ox oy ox oy 

corresponds to the invariant differential equation 

y" = y~^- (4.18) 
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This equation can be obtained from the Lagrangian function 

= / - (4.19) 

y 

which admits all three operators: 

prXiL + LL>(6) = 0; 

prX2L + LD(e2) = 0; (4.20) 

wX^L + LD{i^)^2y'y^D{y-'). 
Consequently, the symmetries yield the following first integrals 



Ji = y" + - = A\ J2 = 2-- 2{y - y'x)y' = 25°, 

y2 y2 



J3^^ + (y-xyr = c' 
y 



(4.21) 



Using the integrals and we write the general solution y{x) as 

AV = (A^'x - By + 1. (4.22) 

We see that the third integral, denoted J3 is not needed, is not useful and 
indeed, is not independent. The integrals Ji, J2 and J3 are related as follows: 

-JiJ3 + l = 0. (4.23) 

Now let us consider the discrete case. Let us consider the discrete Lagrangian 
function 

c^-yl-—. (4.24) 

yy+ 

which admits the same symmetries as the continuous one: 

prXi£ + £D(ei) = 0; 

+h 

prX2>C + £D(6) = 0; (4.25) 
prX3£ + £L»(e3) = L»(y2). 

+h +h 



20 



The Lagrangian generates the invariant quasiextremal equations: 



— : 2{y,^ - ys) = — h „ > 

dy y^y+ y^y_ 

(4.26) 

{y.y + — -{y,)'- — = 0. 

dx yy+ yy. 

The quasiextremal equations have three functionally independent first inte- 
grals 

h = yx^ H = A l2 = + 2y^{x+y^ - y+) = 2B, 

yy+ yy+ 

(4.27) 

^3 = — - + {x+y^ - y+f = C. 

yy+ 

In the discrete case the integrals Ji, I2 and I3 are independent. Eq. ()4.23j) no 
longer holds and instead we have 

I2V . . . I f h 



In order to integrate the system ()4.2(ij] we will use three first integrals A, B 
and the one in eq. ()4.28|) . namely 

e. (4.29) 



yy+ 

Eliminating y^, x+ and ?/_|_, we obtain the solution 

V = {Ax - 5)2 + 1 - ^. (4.30) 

This agrees with the continuous limit ()4.22|) up to order e^. 

Calculating y^ from eq. ()4.30|) and substituting into the expression for A in 
eq. ()4.27|) we obtain a two-point difference equation for x and hence we obtain 
the lattice. In this case the equation has the form of a fractional linear mapping, 
i.e. it is a discrete version of the Riccati equation (with constant coefficients). 

Explicitly we obtain 

aXn -\- j3 f A r>^\ 

Xn+l = -J 4.31 

« = 1-55-ie^ P = ^(1 + B2_i^2). 

^ (4.32) 
7 = -eA; S = l + eB-\e'^. 
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We see that the coefficients in the discrete Riccati equation ()4.31|) satisfy 

a5- (3-i = l, a + (5 = 2 - el (4.33) 

Like the continuous Riccati equation, eq. ()4.31|) can be hnearized. To do this we 
introduce a hnear system 



n 



Un+l \ ^ ( (3\ f U. 

Vn+1 ) \ 6 ) \ Vn 



(4.34) 



If u and V satisfy eq. ()4.34|) . then 



u 



x = - (4.35) 

V 

will satisfy equation ()4.3H) . Eq. ()4.34|) can be solved by standard methods. In- 
deed, both u and v must satisfy 

Un+2 - (« + 5)un+i + {a6 - Pj)Un = 0. (4.36) 

The characteristic equation 

- (a + S)X + {aS - = (4.37) 

is obtained by putting Un = X"'- 

In view of eq. ()4.33j) the roots of eq. ()4.37j) are complex. The final result is 
that the solution of eq. ()4.31|) is 



1 / £2 B 
" A V ~ T A ' ^ ^^^^ 



where p is an integration constant and u satisfies 



tanuj = ^—-^^l-- (4.39) 

Eqs. ()4.30|) and ()4.38p provide an explicit analytic solution of the system 
()4.26|) . It is the general solution and involves 4 constants: A, B, e and p. It 
follows from eq. ()4.38|1 that the independent variable x varies between — and + 
infinity as un + p varies between ±|. 



5 Integration of difference schemes with two variational 
symmetries 

5.1 The method of perturbed Lagrangians 

We have mentioned in Section 3.2 that a two-dimensional group of Lagrangian 
symmetries is always sufficient to reduce the original system of two three-point 
equations to a single three-point equation for the independent variable alone. 
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Here we shall show that in some cases we can do better. Using a different 
approach, we will actually obtain a complete solution of a difference scheme 
approximating a differential equation with a Lagrangian, divergence invariant 
under a two-dimensional symmetry group. 

The case we shall consider is the equation 1)3.211) and hence the two-dimensional 
Abelian group 1)2,1 corresponding to the algebra ()3.20p . We shall make use of 
the fact that the Lagrangian is not unique. Indeed we will consider three differ- 
ent Lagrangians, all having the same continuous limit ()3.22|) . Instead of writing 
the Lagrangian ()3.2fij) in the discrete case, we shall use a family of Lagrangians, 
parametrized by two constants, a and (3: 

C = aGiy,) + Py + (l - P)y+, a ^ I, < P < 1. (5.1) 

Each Lagrangian provides its own quasiextremal system 

a [-G'{y,) + G'{y,)] + (3h+ + (1 - (3)h. = (5.2) 



a [y.G'iy,) - y,G'{y,) - G{y,) + G{y,)]- (3{y -y^) - {l-(3){y+-y) = (5.3) 

We shall view one Lagrangian, with 03 = 1 and jS^ = 0.5 as the basic one, the 
other two as its perturbations. 

Each Lagrangian in the family is divergence invariant under Xi = dx and 
X2 = dy and hence provides two first integrals of the corresponding quasiextremal 
equations ()5.2|) and ()5.3p : 

a [-yxG\y,) + Giy,)] + y + {1 - P)h+yx = A (5.4) 

aG'{yx) -x-ph+ = B. (5.5) 

Let us now choose three different pairs (aj,/?j). They provide six integrals 
(and six quasiextremal equations). We shall show that by appropriately fine 
tuning the constants and jSi and choosing some of the constants Ai and Bi we 
can manufacture a consistent difference system, representing both the equation 
and the lattice. Moreover, we can explicitly integrate the equations in a manner 
that approximates the exact solution obtained in the continuous limit. 

Let us take one equation of the form ()5.4p and two of the form ()5.5p . In these 
three equations we choose cts = 1, P3 = 0.5 and B2 = B3 = B. We then take the 
difference between the two equations involving B to finally obtain the following 
system of three two-point equations: 

«i l-y.G'iy,) + G{yx)] +y+{l- /3i)M. = A (5.6) 

G'iy,) -x-h+ = B (5.7) 
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:i-«2)G"(y,)-(^-/32)/i+ = 0. (5i 



From eq. ()5.7p and ()5.8|) we have 



where we have put 



G'iy.) = ""^^l^^^ (5.9) 



x+-{l + e)x-eB = (5.10) 
2(1 -«2) 



(5.11) 



a2 - 2p2 

The continuous hmit will correspond to e — 0. 

Eq. ()5.10|) coincides with eq. ()4.12|) obtained using three Lagrangian symme- 
tries in a special case. Here it appears in a much more general setting. The 
general solution of eq. ()5.10p 

Xn = {xo + B){1 + e)'' - B (5.12) 

depends on one integration constant xq. This solution gives a lattice satisfying 
h- > and h+ > for xo > —B if e > and for xq < —B if e < 0. For the other 
cases, namely xq < —B if £ > and for xq > —B if £ < 0, formula (j5.12|l gives a 
lattice with a reverse order of points: < and < 0. 

Using ()5.12p and ()5.9|) . we can express i/x in terms of x. We have 

G'(i/J= (l + |)(i? + a;). (5.13) 

Denoting the inverse function of G'{y.j) as iJ, we have 

?/, = i/[(l + |)(S + x)]. (5.14) 

Using ()5.6p and ()5.14p . we can now write the general solution of the system ()5.6|) . 
dSHj) and dEHl) as 

y{x) = A-aiG{H) + {x + B)H, (5.15) 

where we have put 

«i(l + |) -(l-/?i)£ = l. (5.16) 

The value of ai, still figuring in the solution ()5.15|) . must be so chosen as 
to obtain a consistent scheme. Indeed, x„ and given in eq. ()5.12p and ()5.15p 
will satisfy the system (j5.6|) . (j5.7|) and (j5.8|) . We must however assure that 
of eq. (j5.14|) and y^ = {y-n+i — yn)/{.Xn+i — ^n) coincide. A simple computation 
shows that this equality requires that ai should satisfy 
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This equation is consistent only if the right hand side is a constant (independent 
on n). The constants and can depend upon the constant e and for e ^ 
we must have ai, 0:2 ^1; Pi, P2 ^ 0.5. 
From eq. ()5.8j) we have 

This expression must vanish for e 0. To achieve this while respecting eq. ()5.1H] 
we put 

1 

a2 = l + e^, (32 = - + e + -. (5.19) 

Eq. ()5.1ip is satisfied exactly and we have 

h, 2e 



We can view eq. ()5.12|) and ()5.15p as the general solution of the following three 
point difference scheme: 

G'{y.) - G'{y,) - = 0, 

(5.21) 



G'{y,) G'{y,y 

The system fj5.21|) is invariant under the group corresponding to D^,!- Strictly 
speaking, this is not a quasiextremal system, since it can not be derived from any 
single Lagrangian. The arbitrary constants A, B and e come from three first inte- 
grals (|5.6p . (j5.7|) and (j5.2(Jj) that are associated with three different Lagrangians. 

We have not proven that eq. (j5.17p is consistent for arbitrary functions G{yx)- 
We shall however show below that in at least two interesting special cases the 
above integration scheme is consistent. 

The results of this section can be summed up as a theorem. 

Theorem 5.1 The ODE 1)3.211) obtained from the Lagrangian ()3.22|) can be 
approximated by the difference system ()5.21|) . If ai of eq. ()5.17p is constant, 
then the general solution of this system is given by 

X„ = (xo + 5)(l + £)"-S, 

y{xn) = A-aiG{Hr.) + {xn + B)Hn, ^ ' 

where A, B, e and xq are arbitrary constants. For e 0, y{xn) agrees with the 
solution of the ODE (jSm))- 

As applications of this theorem let us consider two different equations, each 
invariant under a three-dimensional group with D2,i as an invariant subgroup. In 
both cases the Lagrangian is only divergence invariant under the subgroup J^2,i- 
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5.2 A polynomial nonlinear ity. 

d d d d 1 

D3,i: ^1 = ^, ^2 = Tr, ^3 = x— + %— , A; ^0,-, ±1,2 (5.23) 
ox ay ox oy 2 

This algebra for k = —1 was treated in Section 4, now we consider the generic 
case. We take 

Giy.) = ^^y.'^-i (5.24) 

and hence 

G'{y,) = {k- l)y,kh = (^1 + (a; + B). (5.25) 
Eq. (|5.14p reduces to 

= Hn{x) = (l + -) (5.26) 



and we have 



Substituting into (j5.17|) . we find 

_ t(l +,)((!+ ef-'-l) 
- (*-l)(l + i)((l 

so that we have ai = 1 + 0(£^). 

Thus, ai is a constant, close to ai = 1 for e <^ 1. The solution yn of ()5.22p 
specializes to 

+ (/t- 1)^-1 (l+e)^-l' ^ ^ ^ 

This agrees with the solution (lOl) of the ODE up to 0(e2). 

It is interesting to note that for k = —1 ai becomes independent on e and we 
obtain ai = 1, /3i = 0.5. The solution ()5.29p provides us with the solution ()4.16p . 
which was obtained in Section 4 with the help of a different method. 



5.3 An exponential nonlinearity 

We consider another three-dimensional group and its Lie algebra, namely: 

d d d d 

03,3: ^i = 7r' ^2 = ^, X3 = x— + {x + y)—. (5.30) 
ox oy ox oy 

The corresponding invariant ODE is 

y" = exp{-y') (5.31) 
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and can be obtained from the Lagrangian 

L = exp(y') + y. 

We have 

prXiL + LD{^i) = 0; 
prX2L + LZ}(6) = l = D{xy, 

The corresponding first integrals of eq. ()5.31|1 are 

eyipiy')i^ - y') + y = A exp{y')-x = B. 
Finally, the general solution of eq. is 

y={x + B){ln{x + B) - 1) + A. 



(5.32) 
(5.33) 

(5.34) 
(5.35) 



Now let us consider the discrete case, following the method of Section 5.1. 
We have 

Giyx) = exp(y^) 

and hence 

G'iy^) = exp(y^) = (x„ + B) {l + |) 
Hn = y, = ln{xn + 5) + In (1 + I) 

Substituting into eq. ()5.17|) . we find 



(5.36) 
(5.37) 



ai 



[l+e) ln(l +£) 
e{l + l) 



(5.38) 



so that «! is indeed a constant and moreover we have ai = 1 + O^e"^). 
The solution y{x) on the lattice given in eq. (j5.22|) is 



yn = A+{xn + B) \n{xn + B) + {xn + B) 



£^ (l + e)ln(l + £) 



(5.39) 



This agrees with the solution ()5.35|1 of the ODE ()5.31|1 up to 0{e^ 



6 Concluding remarks 

We see that variational symmetries, and the first integrals they provide, play a 
crucial role in the study of exact solutions of invariant difference schemes. Much 
more so than in the theory of ordinary differential equations. 

The procedure that we followed in this article can be reformulated as follows. 
We start from the continuous case where we know a Lagrangian density L{x, y, y'), 
invariant under a group Gq of local point transformations, i.e. satisfying condition 
(j2.6|) . or (j2.7|) . We hence also know the corresponding Euler-Lagrange equation, 
invariant under the same group, or a larger group containing Gq as a subgroup. 



27 



Let us assume that we can approximate this Lagrangian by a "discrete La- 
grangian density" L{x,y,x+,y+) invariant under the same group Gq. Even in 
the absence of any symmetry group, the Lagrangian will provide us with the 
quasiextremal equations ()2.2H) . i.e. with a discrete Euler-Lagrange system. This 
system can be identified with the difference system p. HI . p.2|l . 

If the Lagrangian is invariant under a one- dimensional symmetry group, we 
can reduce the quasiextremal system to a three-point relation for x alone, plus 
a "discrete quadrature" for y (see Section 3.1). If the symmetry group of the 
Lagrangian is two-dimensional, we can always reduce the quasiextremal system 
to one three-point equation for x alone, and write the solution yn{x) directly (see 
Section 3.2). 

If the invariance group of the Lagrangian is (at least) three-dimensional then 
we can integrate the system explicitly (Section 4). 

Finally, we have shown that if the symmetry group of the Lagrangian is two- 
dimensional, but the quasiextremal system has a third (non-Lagrangian) symme- 
try, we can also integrate explicitly. 
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